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FOREWORD 


This  report  was  prepared  by  Dr.  T.  Y.  Yang,  Visiting  Scientist 
In  the  Optimization  Group,  Analysis  and  Optimization  Branch,  Structural 
Mechanics  Division,  Air  Force  Flight  Dynamics  Laboratory  (AFFDL/FBR), 

Air  Force  Wright  Aeronautical  Laboratories,  Air  Force  Systems  Command, 
Wright-Patterson  Air  Force  Base,  Ohio.  This  work  was  performed  under 
Project  Nr.  2307,  "Research  In  Flight  Vehicle  Dynamics",  Task  Nr.  23070S, 
"Basic  Research  In  Structures  and  Dynamics". 

This  work  was  also  performed  under  Contract  F33615-76-C- 31^5  by 
T.  Y.  Yang,  Professor  of  School  of  Aeronautics  and  Astronautics,  Purdue 
University,  West  Lafayette,  Indiana  47907. 

The  technical  report  was  released  by  the  Author  In  October  1977. 

The  report  covers  work  conducted  from  June  1976  through  October  1977. 
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SECTION  I 


INTRODUCTION 

Since  van  der  Neut  (Reference  1)  demonstrated  the  influence  of  eccentricity 
of  stiffeners  on  the  buckling  strength  of  stiffened  cylindrical  shell  in  1947, 
numerous  papers  on  the  buckling  analysis  of  stiffened  cylindrical  shells  using 
the  smeared-out  technique  were  published  (References  2-10).  Such  studies  were 
motivated  by  the  structural  efficiency  and  practical  advantages  of  the  wall- 
stiffening.  Reference  3  showed  that  when  the  smeared-out  method  was  used  for 
cylinders  with  large  number  of  stiffeners,  good  correlation  between  experi¬ 
mental  and  theoretical  results  in  buckling  load  was  obtained.  However,  this 
may  not  be  the  case  for  cylinders  with  small  number  of  stiffeners. 

MacNeal  ^  al_.  (Reference  11)  treated  ring-stiffened  cylindrical  shells  by 
considering  rings  as  discrete  elements.  Wang  and  Lin  (Reference  12)  and  Singer 
and  Haftka  (Reference  13)  analyzed  the  buckling  of  cylindrical  shells  with  dis¬ 
crete  stringers  and  showed  the  differences  in  results  between  the  smeared-out 
and  the  discrete  methods. 

Egle  and  Sewall  (Reference  14)  derived  the  frequency  equations  for  the 
cylindrical  shells  with  discrete  stringers  and  rings.  Useful  information  on 
circumferential  modes,  which  are  in  common  for  both  vibration  and  buckling  prob¬ 
lems,  was  given.  Their  numerical  results  were  limited  to  the  study  of  the  ef¬ 
fect  of  discrete  stringers  on  natural  frequencies.  McDonald  (Reference  15)  solved 
a  problem  on  the  free  vibration  of  stiffened  cylindrical  shells  by  considering 
rings  as  discrete  elements. 

Wang  (Reference  16)  performed  the  deformation  and  stress  calculations  for 
orthogonally  stiffened  cylinders  under  internal  pressure  by  treating  the  stif¬ 
feners  as  separate  elements.  Soong  (Reference  17)  analyzed  the  general  instability 
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of  orthotropic  stiffened  cylinders  with  intermittently  attached  stiffeners  by 
considering  the  discreteness  of  the  stiffeners.  The  buckling  patterns  assumed 
were  symmetric  circumferential  and  symmetric  longitudinal  modes.  The  results 
calculated  on  the  basis  of  those  displacement  functions  for  the  cylindrical 
shells  with  many  stiffeners  agreed  well  with  experimental  data  and  smeared-out 
calculations . 

In  this  report,  the  discreteness  of  the  stiffeners  is  considered  in  the 
buckling  analysis  of  orthogonally  stiffened  cylinders  with  classical  simply 
supported  edge  conditions.  It  is  shown  that  when  all  the  stringers  are  identical 
and  equally  spaced,  the  buckling  equations  can  be  reduced  to  two  sets  of  equations, 
one  set  yields  the  antisymmetric  circumferential  modes  and  the  other  yields  the 
symmetric  circumferential  modes.  Each  set  can  then  be  uncoupled  into  several 
circumferentially  related  equations.  It  is  also  shown  that  when  all  rings  are 
identical  and  equally  spaced,  the  buckling  equations  can  be  reduced  to  two  sets  of 
equations,  one  set  yields  the  antisymmetric  longitudinal  modes  and  the  other 
yields  the  syi^imetric  longitudinal  modes.  The  uncoupled  submatrices  are  further 
reduced  by  partitioning  and  substitution.  Effort  is  made  to  preserve  the 
sparseness  of  the  matrices  in  order  to  save  computing  time  and  storage  by  using 
a  special  compact  storage  scheme  (Reference  20). 

The  large  sets  of  general  eigenvalue  buckling  equations  are  solved  by  the 
Ritz  iteration  method  combined  with  Chebyshev  procedure  (Reference  18).  The 
accuracy  and  efficiency  of  the  method  are  evaluated  through  an  example. 

Examples  of  buckling  predictions  for  discretely  stiffened  cylinders  are 
performed  and  results  are  compared  with  available  analytical  and  experimental 
solutions  (Reference  19). 

The  problems  of  minimum  weight  design  of  stiffened  cylindrical  shell  under 
buckling  constraints  have  long  been  of  interest  to  the  aerospace  structural 
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designers.  The  studies  of  such  problems  have  been  extensive  (References  21-37). 

Before  1967,  the  studies  were  based  on  the  assumption  of  simultaneous  failure 
modes  and  the  optimization  was  achieved  by  parametric  studies  (References  21-23). 
The  structural  syntheses  for  the  subject  problems  were  performed,  among  others, 
by  Kicher  (Reference  24),  Schmit,  Morrow,  and  Kicher  (Reference  25),  and  Pappas  and 
Amba-Rao  (Reference  28).  Kicher  treated  the  problem  by  using  the  constrained 
gradient  method.  Schmit,  Morrow,  and  Kicher  applied  a  Fiacco-McCormick-type 
penalty  function  formulation  to  transform  the  basic  inequality  constrained  min¬ 
imization  problem  into  a  sequence  of  unconstrained  minimization  problem.  Pappas 
and  Amba-Rao  used  a  direct  search  algorithm  with  an  interior-exterior  penalty 
function  formulation. 

Thornton  performed  the  synthesis  of  stiffened  conical  shells  using  the  ex¬ 
terior  penalty  function  method  with  least  square  approximation  (Reference  32). 

Kunoo  and  Yang  carried  out  the  minimum  weight  design  of  cylindrical  shells  with 
orthogonal  stiffeners  subjected  to  uniform  axial  compressive  or  bending  load 
by  the  method  of  steepest  descent  (Reference  37). 

Hofmeister  and  Felton  showed  that  the  waffle  plate  with  multiple  rib  sizes 
compared  more  favorably  to  sandwich  plates  than  those  with  a  single  stiffener 
size  (References  38-39). 

In  all  the  references  (References  21-37)  concerned  with  the  weight  opti¬ 
mization  of  the  stiffened  shells,  it  appears  that  all  the  buckling  predictions 
were  based  on  the  formulations  using  smeared-out  technique  to  treat  the  stiffen¬ 
ers.  The  main  reason  for  not  using  discrete  design  and  discrete  formulation  is 
that  they  result  in  large  sets  of  eigenvalue  equations  that  are  extremely  ex¬ 
pensive  to  solve  repeatedly  during  optimization. 

In  this  report,  the  exact  buckling  formulations  for  the  cylindrical  shells 
with  discrete  stiffeners  are  derived.  By  assuming  equally  spaced  and  identical 
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discrete  stringers  and  rings,  the  buckling  formulations  can  be  uncoupled  into 
simpler  ones  based  on  their  associated  mode  shapes.  Then  following  some  matrix 
partitioning  and  reduction,  some  efficient  use  of  matrix  sparsities,  and  the  use 
of  efficient  eigenvalue  solution  procedures,  the  huge  sets  of  eigenvalue  equations 
can  be  reduced  and  become  manageable.  It  is  demonstrated  in  a  typical  example 
of  cylinder  with  discrete  stiffeners  that  approximately  400  seconds  of  central 
processing  time  are  required  for  computing  one  lowest  eigenvalue  by  a  CDC  6500 
computer.  Such  problem  was  governed  by  750  eigenvalue  equations  in  the  original 
coupled  formulation.  The  computing  time  of  400  seconds  is,  however,  still  formidably 
long  in  the  minimum  weight  design  where  repeated  computations  of  the  lowest 
eigenvalues  are  required. 

To  circumvent  this  difficulty,  it  is  suggested  here  to  assume  extremely 
simple  displacement  functions  in  the  buckling  formulation.  With  such  simplicity, 
the  computing  time  required  for  repeated  eigenvalue  computations  in  the  optimiza¬ 
tion  process  becomes  small.  To  ensure  that  the  optimized  design  is  feasible, 
the  exact  buckling  load  based  on  the  exact  discrete  formulation  is  finally  com¬ 
puted  using  the  optimized  design  variables  to  check  the  buckling  load  obtained 
by  the  approximate  method. 

A  method  of  multipliers  recently  developed  by  Schuldt  (Reference  40)  for 
mathematical  programming  problems  is  extended  to  the  present  structural  problems 
for  efficient  minimum  weight  design. 

Two  design  examples,  one  with  the  smeared-out  and  one  with  both  the  smeared- 
out  and  the  discrete  stiffeners,  are  performed  to  illustrate  the  present  simplified 
buckling  analysis  method  and  the  method  of  multipliers  for  optimization. 
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SECTION  II 


BUCKLING  FORMULATIONS 


1.  BASIC  ASSUMPTIONS 

For  the  shell  skin,  it  is  assumed  that  the  Kirchhoff  and  Love  assumptions 
for  thin  elastic  shell  hold. 

For  the  stiffeners,  the  following  assumptions  are  made. 

(a)  The  shell  is  integrally  stiffened  by  longitudinal  stringers  and  cir¬ 
cumferential  rings  for  both  smeared-out  and  discrete  types. 

(b)  The  smeared-out  stiffeners  are  thin  and  deep  such  that  the  uniaxial 
stress  condition  prevails  and  the  torsional  stiffness  can  be  neglected. 

2.  TOTAL  POTENTIAL  ENERGY  EXPRESSIONS 

The  total  potential  energy  for  the  system  as  shown  in  Figure  1  is  defined  as 

TT  =  U^  +  U2  +  U3  +  U^  -  -  W2  -  W3  (1) 

where  U-j  is  the  strain  energy  for  the  shell  skin;  U2  is  the  strain  energy  for  the 
smeared-out  orthogonal  stiffeners;  U3  and  are  the  strain  energies  for  the 
discrete  stringers  and  rings,  respectively;  W.j,  ^2*  and  W3  are  the  potential  energies 
due  to  the  middle-surface  load  for  the  shell  with  the  smeared-out  stiffeners, 
the  discrete  stringers,  and  the  discrete  rings,  respectively. 

The  strain  energy  expressions  written  in  terms  of  strains  are  available 
collectively  in,  for  example,  References  2,  10,  and  14.  These  expressions  can  be 
derived  in  terms  of  displacements  by  using  the  following  strain-displacement 


relations 


^2  -  .0/.  w/r  -  ^vv,0  -  w,00;/r 

=12  '  ’.x  *  V"  •  ^"‘“ixe  * 

It  is  noted  that  the  terms  zv  ./r^  and  zv  /r  are  neglected  if  Donnell  type 

|0  >  A 

theory  is  used. 


The  potential  energy  versus  displacement  expressions  take  the  following 


forms , 
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3.  DISPLACEMENT  FUNCTIONS 

The  displacement  functions  which  satisfy  the  classical  simply-supported 
boundary  conditions  that  w  »  -  N^^  *  v  «  0  at  x  *  0  and  L  may  be  represented 


'  u"  •  k'-l 
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n=i 


mirx 


4.  EQUILIBRIUM  EQUATIONS  FOR  BUCKLING  ANALYSIS 

Substituting  the  displacement  functions  (6)  into  the  expressions  for  the 
total  potential  energy  and  then  using  the  minimum  potential  energy  principle, 
the  equations  of  equilibrium  are  obtained  as 
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where  the  coefficients  in  the  two  symmetric  matrices  are  defined  as 
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As  to  the  double  signs,  the  upper  and  lower  ones  are  for  the  inner  and  outer  stiffened 
shells,  respectively.  The  Kronecker  delta  function  is  defined  as  6^^  «  0  if  i  j 
and  *  1  if  i  *  j.  For  the  matrix  elements  with  subscripts  1  and  2,  it  is  defined 

that 

m  •  m^  »  '"j  *  ^ 

n  »  n^  «  i  -  [i/(k^  +  l)]i  '’j  “  ■*■ 

I 

For  the  rest  of  the  matrix  elements,  Eq.  (9)  also  applies  except  that  k^  is  replaced 
in 

I  • 

The  order  of  the  matrices  in  Eqs.  (7)  is  (4k  k_  +  2k  k).  The  submatrices 

m  n  m  n 

C  and  D  are  functions  of  the  properties  of  the  discrete  stringers  only.  If 
there  is  no  stringers,  Eqs.  (7)  reduce  to  two  sets  of  equations  which  result  in  the 
same  buckling  load  but  two  circumferential  modes  with  90**  phase  difference. 

The  uniform  axial  load  and  uniform  external  pressure  are  considered  here. 

They  are  related  to  the  internal  stresses  by  using  the  assumption  of  either  uniform 
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stresses  or  uniform  strains  in  the  system. 

The  stresses  based  on  the  uniform  stress  assumption  are 

%  *  ^s  °e  “  ^  ^r  ^  ^2^  ^10) 

The  stresses  based  on  the  uniform  strain  assumption  are 

°x  ”  ^^^^2  +  “r 
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where  H  -  (H^  +  +  1)(H2  +  +  1)  -  v^(H^  +  H^)(H2  +  H^);  -  h^^/h;  H2  «  h^g/h; 

Hs  =  E^h^^/Eh;  and  *  E^h^^/Eh. 

The  axial  load  and  external  pressure  have  the  relation  that  >»  pr/2. 

5.  EQUALLY  SPACED  AND  IDENTICAL  DISCRETE  STRINGERS  AND  RINGS 

If  the  stringers  are  identical  in  size,  the  matrix  elements  with  triple  primes 
in  Eqs.  (7)  include  a  common  term, 

L  *^5 

"l"j  *  kh  '’i^k  "j®k  ^’2) 

For  equally  spaced  stringers,  =  2Trk/n^  and  Eq.  (12)  becomes 
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Thus  It  Is  shown  that  for  the  case  of  Identical  and  equally  spaced  stringers,  the 
submatrix  elements  with  triple  primes  In  E,s.  (7)  uanish.  Egs.  (7)  are  thus  un¬ 
coupled  Into  two  sets  of  equations.  One  set  Is  related  to  the  amplitudes  u'.  v'. 
and  W  and  yields  the  antisymmetric  circumferential  buckling  modes.  The  other  Is 

related  to  the  amplitudes  u".  v".  and  w"  and  yields  the  symmetric  circumferential 
buckling  modes. 

The  elements  In  the  submatrices  with  single  and  double  primes  contain  two 
common  terms 
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using  the  similar  procedures  as  used  In  deriving  Eqs.  (13  and  14).  Eq.  (15) 
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For  Identical  and  equally  spaced  rings,  the  buckling  equations  (7)  Include 
the  following  terms, 

with  X,^  -  kL/(n^  +  1)  (18) 

Similar  derivations  give 
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1 

Further  examination  of  Eqs.  (20)  to  (23)  reveals  that  If  m^-m^  0  (mod  2)  and 

m.+m.  ^  0  (mod  2), 

*  J 


This  simple  conclusion  leads  to  the  uncoupling  of  the  buckling  equations  (7) 
into  two  sets,  one  gives  the  symmetric  longitudinal  buckling  modes  and  the  other 
gives  the  antisymmetric  longitudinal  buckling  modes. 


The  equations  derived  here  Involve  no  summations,  the  computation  Is  thus 


greatly  simplified. 


6.  FURTHER  CONSIDERATION  OF  IDENTICAL  AND  EQUALLY  SPACED  DISCRETE  STRINGERS 
AND  RINGS 

A  typical  stiffness  matrix  element  In  the  buckling  equations  (7)  may  be 
written  In  the  symbolic  form 

c(1,j)  «  6..A  +  „  (B  „  +  D  ^  (E  Q„  „  +  F  R,  „  ) 

ij  m^mj  n^n^  n-n^  n^n^  m^m^  'in^mj 


Eq.  (25)  Indicates  that  odd  numbers  of  m^  and  the  even  numbers  of  m^  are  not 
related.  In  other  words,  If  a  row  In  the  buckling  matrix  corresponding  to  some 
odd  number  of  m^ ,  then  all  the  elements  In  that  row  that  correspond  to  the  even 
numbers  of  mj  are  zero.  For  a  sample  case  with  two  stringers  (n^  *  2),  the  non¬ 
zero  elements  In  the  matrix  are  shown  In  Fig.  2a.  By  properly  Interchanging  rows 
and  columns  In  Fig.  2a,  the  uncoupled  version  of  matrix  Is  obtained  and  shown  In 
Fig.  2b.  Fig.  2b  shows  that  the  submatrix  that  corresponds  to  the  odd  numbers  of 
m^  and  mj  results  In  a  symmetric  longitudinal  buckling  mode,  whereas  the  other 
submatrix  that  corresponds  to  the  even  numbers  of  m^  and  mj  results  In  an  antisymmetric 
longitudinal  buckling  mode. 

Extending  this  Idea  to  the  circumferential  buckling  mode,  Eqs.  (16)  can  be 


reduced  as  follows 


M  *  N 
"l^j  "l"j 
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n.  +  n.  ^  0  (mod  n  ) 

1  J  5 


Eqs.  (26)  Indicate  that  If  some  row  In  the  buckling  matrix  corresponds  to  n^,  then 
all  the  elements  In  that  row  that  corresponds  to  the  number  n.  which  satisfies 

J 

n^  -  nj  ^  0  (mod  n^)  and  n^  +  nj  |  0  (mod  n^)  are  zero.  In  other  words,  the 
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(a)  The  nonzero  elements  In  a  typical  buckling  matrlxt 

(b)  The  two  uncoupled  submatrices  with  symmetric  and 
antisymmetric  longitudinal  modes;  and 

(c)  The  four  further  uncoupled  submatrices  with 
circumferentially  related  mades. 


elements  with  and  nj  that  satisfy  n.  -  n^  =  0  (mod  n^)  and  n^  +  nj  =  0 

(mod  n  )  are  not  zero.  For  example,  if  n  =  8,  the  elements  that  correspond 

S  5 

to  the  combinations  with  (n^.nj  =  1,  7,  9,  15,  17,  ....),  (n^,nj  »  2,  6,  10,  14, 
....),  (n.,n.  =  3,  5,  11,  13,  19,  21,  ....),  (ni,n.  *  4,  12,  20,  28,  and 

I  J  •  J 

(n.,n.  =  0,  8,  16,  24,  32,  _ )  have  non-zero  values.  These  relations  may  be 

*  0 

illustrated  in  Figs.  2a,  2b,  and  2c  for  the  case  of  n^  =  2. 

These  procedures  can  be  summarized  as  follows.  The  whole  set  of  buckling 
equations  is  first  uncoupled  into  two  sets:  one  related  to  the  symmetric  circum¬ 
ferential  modes  and  the  other  related  to  the  antisymmetric  circumferential  modes. 
Each  of  these  two  sets  is  then  uncoupled  into  two  sub-sets:  one  related  to  the 
symmetric  longitudinal  modes  and  the  other  related  to  the  antisymmetric  longitudinal 
modes.  The  lowest  eigenvalue  that  comes  from  these  four  sub-sets  is  the  buckling 
load  of  the  whole  system. 


SECTION  III 


NUMERICAL  METHOD  FOR  BUCKLING  EQUATIONS 


As  shown  in  Eq.  (14),  for  the  case  of  equally  spaced  and  identical  stringers, 
the  buckling  equations  (7)  can  be  reduced  into  two  sets;  one  related  to  the  anti¬ 
symmetric  circumferential  modes  and  the  other  related  to  the  symmetric  circumfer¬ 
ential  modes.  Both  sets,  although  different  in  content,  can  be  described  in  the 
same  symbolic  form. 


("V  'X,  ^ 

/  'X/  \ 

0  0  0 

0 

X  ^  ^ 

'V/ 

% 

^23 

- 

0  0  0 

< 

V 

-  B  1 

0 

^  a;  “V 

'U 

^33^ 

“  “33  J 

1  wj 

0 

(27) 


In  order  to  save  computer  storage  and  computing  time  by  a  special  compact 
storage  scheme,  effort  to  preserve  the  original  form  of  the  sparseness  of  the 
buckling  submatrices  is  made. 


Through  some  manipulations,  Eqs.  (27)  are  reduced  to 

Cl  =  A  Cl 

where  A  =  ^33  - 

6  =  A  a 

a  =  a  -  a  a 

c  1*22 


23 


'23 


(28) 


(29) 
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For  uniform  axial  compression  and  external  pressure,  matrices  A  and  A  are  both 
positive  definite  and  symmetric.  Because  only  the  smallest  eigenvalue  A  is  of 
interest,  the  Power  method  may  be  used.  Unfortunately,  for  the  present  class  of 
problems  the  eigenvalues  are  poorly  separated  and  the  convergence  is  slow. 


An  effective  algorithm  called  Ritz  iteration  was  developed  by  Rutihauser  (Refer¬ 
ence  18)  for  standard  eigenvalue  problem.  It  can  be  extended  to  solve  the  eigenvalue 
equations  (28)  via  the  following  procedure. 

Matrix  A  is  first  decomposed  as  the  product  of  a  lower  triangular  matrix  t 
and  its  transpose  by  means  of  Cholesky  factorization.  Eq.  (28)  can  then  be 

derived  as 


)-‘i  =  i  it  with  it .  t^ii 


(30) 


Matrix  Sf  is  the  initial  matrix  in  iteration,  whose  columns  are  assumed  as  or¬ 
thonormal  and  the  order  of  the  matrix  is  pxq  with  p  being  the  order  of  matrix 
A  and  q<p.  The  iterative  procedure  is  described  by  a  typical  cycle  from  step 


k-1  to  k. 

(a) 

4 

'  ^  -  (t^)-^Vi 

1 

II 

u 

o 

(31) 

or  t  \  -  !) 

(b) 

(32) 

(c) 

Solve  for  the  eigenvalues  of  the  qxq  matrix 

(d) 

\  ■  Wk“ 

(33) 

where  is 

the  matrix  of  eigenvectors  of  matrix  and  is  the  diagonal 

matrix 

whose  elements  are  the  square  roots  of  the  eigenvalues  of  matrix  in  decreasing 
order.  The  eigenvalues  of  Eq.  (28)  are  calculated  by  the  Rayleigh  quotient 
formula 


X  =  1  t 


xA  ^B(J)  ^x. 


(34) 
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(35) 


and  the  associated  eigenvectors  are  given  by  Eq.  (30)  as 

Cf  = 

The  convergence  rate  of  this  method  Is  considerably  Improved  over  the  Power 
Method.  It  Is,  however,  still  slow  for  the  present  problems  because  quite  a  few 
eigenvalues  are  very  close  to  the  minimum  eigenvalue.  Rutlshauser  also  suggested 
that  If  Chebyshev  Iteration  Is  used  In  addition  to  the  Iterative  step  (a)  In  Eq. 
(31)  considerable  Improvement  In  convergence  rate  can  be  achieved,  especially  when 
convergence  Is  slow. 


The  Chebyshev  Iteration  Is  expressed  as 

^h+1  '  e  ^  ^  ^Vl-i  ■  ^\+1-i  ■  Vl-2 


( 1~2 ,3 , . . . , j -1 ) 


(36) 

(37) 


where  e  meets  the  condition  that 


:  2e  ; 


(38) 


The  entire  procedure  applied  here  may  be  described  In  the  following  nine 
steps. 

(1)  Choose  a  set  of  values  In  matrix  5^p(pxq). 

(2)  Orthonormal Ize  the  columns  In  matrix  51^  by  the  Grani-Schmidt  orthogonal Izatl on 
process (Reference  43) . 

(3)  Start  with  h  =  0  and  j  =  2. 

(4)  Compute  j  -  1  steps  of  values  by  using  Eqs.  (36)  and  (37). 

(5)  Orthonormal  Ize  the  columns  of  .  .  . 

h+j-i 

(6)  Perform  one  step  of  Iteration  as  described  In  (a),  (b),  (c)  and  (d)  with 
k  =  h  +  j. 
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(7)  Determine  j  such  that  It  satisfies  the  condition 


j  -  1  <  Arccosh  (10)/Arccosh  {(a^^  -  e)/e}  (39) 

(8)  Check  tolerance  In  the  lowest  eigenvalue  between  two  current  steps. 

(9)  Return  to  step  (4). 

This  procedure  was  successfully  applied  to  the  present  large  eigenvalue 
equations  which  contain  poorly  separated  eigenvalues.  This  method  can  also  cal¬ 
culate  the  equal  eigenvalues  which  may  appear  In  the  vibration  problem. 

The  convergence  characteristics  of  the  Power  Method,  the  Ritz  Iteration  Method, 
and  the  Ritz  Iteration  Method  combined  with  Chebyshev  procedure  were  studied 
through  a  numerical  example.  The  example  was  derived  from  Eqs.  (27)  with 
km  =  kfj^  »  4  and  kn  *  =  8.  The  matrix  given  by  Eqs.  (27)  had  an  order  of 

96  which  was  reduced  to  32  by  using  Eqs.  (28).  The  results  are  summarized  In 
Table  1.  The  Ritz  iteration  method  with  Chebyshev  procedure  appears  to  be  most 
efficient  among  the  methods  compared.  This  method  was  used  in  the  example  studies. 


Table  1.  Comparison  of  Convergence 


Method  1 

Mode  Number 

CDC  6500 

CP  Time 
(sec.) 

Iteration 

Lowest 

2nd 

3rd 

4th 

Exact 

1511.627 

1544.382 

1590.751 

1772.99 

Power 

1511.630 

12.0 

245 

Ritz  Iteration 

1511.629 

1544.393 

1590.759 

1797.46 

7.3 

30  X  4 

Ritz  +  Chebyshev 

1511.627 

1544.382 

1590.751 

1778.94 

4.5 

22  X  4 
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SECTION  IV 


NUMERICAL  RESULTS  FOR  BUCKLING  ANALYSIS 

1.  SIMPLY-SUPPORTED  CYLINDER  WITH  MANY  STIFFENERS  UNDER  UNIFORM  AXIAL 
COMPRESSION  -  BY  BOTH  SMEARED-OUT  AND  DISCRETE  TECHNIQUES 

The  first  example  was  defined  by  the  dimensions  and  properties  that  E  ■  1.06  x 
10^  psi;  V  *  1/3;  a^  ■  a2  •  3  In.;  d^  *  d2  “  0.2  In.;  *  t2  ■  0.125  In.;  h  «  0.05  In.; 
r  »  48  In.;  L  *  50  In.,  and  weight  »  100  pounds.  The  results  are  given  In  Table  2. 


Table  2.  Uniform  Axial  Buckling  Load  P*  for  Example  11) 


Method 

Px 

(kips) 

P-/2irr 

(Ibs/ln.) 

°x 

(psT) 

Mode  Numbers 

m 

n 

(1) 

D1  screte 

808.0 

m 

6 

17 

(2) 

Discrete 

801.3 

6 

17 

(3) 

Smeared-out 

791.5 

WBM 

6 

17 

In  the  discrete  calculations,  16  rings  and  100  stringers  were  used.  The  first 
discrete  calculation  was  based  on  the  displacement  functions  with  600  terms  (k^^^  ■ 

I  I 

k^  *  *  ^'n  *  antisymmetric  circumferential  mode,  and  uniform  strain 

distribution.  The  second  discrete  calculation  was  based  on  five  antisymmetric 
circumferentially  related  modes  and  24  symmetric  longitudinal  modes  (k^  ■  5  with 
n  ■  17,  83,  117,  183,  217;  ■  24  with  m  ■  1,  3,  5,  ...47)  and  uniform  strain 

distribution.  The  smeared-out  calculation  was  based  on  uniform  stress  assumption 
and  the  result  agreed  with  that  given  In  Reference  10.  Table  2  serves  to  verify 
the  correctness  of  the  present  discrete  formulations  and  solution  methods.  Be¬ 
cause  the  torsional  stiffnesses  of  the  stiffeners  were  neglected  In  the  smeared- 
out  calculation,  the  resulted  buckling  load  should  be  slightly  smaller  than  those 
obtained  by  the  discrete  calculations. 


When  the  numbers  of  rings  and  stringers  are  large,  the  buckling  loads  obtained 
by  using  uniform  strain  and  uniform  stress  assumptions  should  be  close.  The  first 
discrete  computation  was  also  performed  by  using  the  uniform  stress  assumption,  a 
buckling  load  of  245.7  kips  v/as  obtained.  Since  the  uniform  strain  assumption 
Introduces  circumferential  tensions  In  the  rings  which  produce  compressive  stresses 
In  the  skin.  It  always  yields  lower  buckling  load  than  the  uniform  stress  assumption. 

If  the  stiffeners  in  this  example  are  neglected,  the  buckling  stress  o  for 

A 

the  monocoque  cylinder  becomes  6750  psi.  The  tremendous  strengthening  effect  due  to 
a  large  number  of  stiffeners  Is  seen. 

2.  SIMPLY-SUPPORTED  CYLINDER  WITH  ONE  RING  AND  FOUR  STRINGERS  UNDER  UNIFORM 
AXIAL  COMPRESSION 

The  geometry  and  the  material  parameters  of  the  second  examples  are  defined 
as:  V  =  1/3;  h  =  0.03  In.;  r  =  48  In.;  L  =  50  In.;  *  d^  =  1.3196  In.;  t^  * 
t^  =  0.8248  In.;  Z  »  1637;  and  weight  *  100  pounds.  The  buckling  stresses  computed 
by  using  the  uniform  stress  assumption  and  either  symmetric  or  antisymmetric  circum¬ 
ferential  mode  are  given  In  Table  3. 

The  converged  buckling  stress  Is  4070  psi.  Table  3  reveals  that  one  should  not 
be  prejudiced  against  the  u  and  v  functions  by  using  more  terms  In  the  w  function. 

The  buckling  calculation  was  then  performed  by  assuming  uniform  stress  and 
either  symmetric  or  antisymmetric  circumferential  mode,  elgher  symmetric  (m  =  1 , 

3,  ...  2k|^  +  1)  or  antisymmetric  (m  =  2,  4,  ...  2k|^)  longitudinal  modes,  and 
circumferentially  related  modes.  The  results  for  buckling  stresses  and  mode  shapes 
are  shown  in  Figs.  3a,  3b,  4a,  and  4b,  respectively.  The  mode  shapes  were  plotted 

by  the  Gould  4800  electrostatic  printer.  The  four  mode  shapes  In  the  figures  are  quite 


Table  3.  The  Buckling  Stresses  For  The  Discretely  Stiffened  Cylinder 
Antisymmetric  circumferential  mode 


CDC-i 


V 

No.  of 
terms 

Buckling  Stress 
0^  {lbs/in2) 

Relative 

Error 

12 

2 

12 

2 

72 

10107.0 

3.6x10"® 

10x4 

4 

14 

2 

14 

2 

84 

8165.6 

8.2x10"® 

12x4 

5 

13 

3 

13 

3  * 

117 

9498.3 

1.8x10"^ 

22x4 

13 

12 

4 

12 

4 

144 

10107.0 

6.5x10"^ 

19x4 

19 

24 

2 

24 

2 

144 

5581.3 

3.4x10"® 

33x4 

25 

26 

2 

26 

2 

156 

5576.0 

1.6x10"^ 

50x4 

36 

28 

2 

28 

2 

168 

5575.0 

6.3x10"^ 

51x4 

43 

30 

2 

30 

2 

180  1 

1 

5574.6 

3.4x10"^ 

58x4 

54 

24 

4 

24 

4 

288 

5581.3 

9.9x10"^ 

79x4 

288 

10 

10 

10 

10 

300 

11319.0 

5.3x10"^^ 

38x4 

142 

10 

11 

10 

11 

330 

10916.0 

1.0x10"^ 

30x4  1 

- 

11 

10 

11 

10 

330 

9910.3 

4.2x10"® 

27x4 

156 

28 

8 

14 

4 

336 

4320.3 

4.7x10"^ 

115x5 

922 

9 

13 

9 

13 

351 

11798.0 

6.5x10"^ 

52x4 

239 

25 

5 

25 

5 

375 

4070.3 

2.0x10"® 

43x4 

234 

26 

6 

26 

6 

468 

4070.3 

4.5x10"® 

59x4 

536 

28 

6 

28 

6 

504 

4070.0 

9.9x10"^  1 

84x4 

632 

26 

7 

26 

7 

546 

4071.6 

3.7x10"®  1 

40x4 

582 

27 

7 

27 

7 

567 

4070.3 

7.9x10"^  I 

114x4 

964 

26 

8 

26 

8 

624 

4070.6 

2.0x10"®  1 

35x4 

901 

Symmetric  circumferential  mode 


12 

2 

12 

2 

71 

17295.0 

immii 

5 

14 

2 

14 

2 

83 

17295.0 

6 

13 

3 

13 

3 

116 

16808.0 

2.3x10"' 

25x4 

14 

12 

4 

12 

4 

143 

11073.0 

2.2x10"^ 

8x4 

15 

24 

4 

24 

4 

287 

5582.3 

4.2x10"^ 

33x4 

114 

25 

5 

25 

5 

373 

4070.3 

3.9x10"^ 

39x4 

271 

10 

17 

10 

17 

539 

8383.0 

1.9x10"^ 

27x4 

457 

Figure  3  Buckling  Mode  Shapes  for  Anti  symmetric  Circumferential  Mode  and 

(a)  Mixed  Longitudinal  Modes  (a  =  4070  psi ;  k  =  28;  k  ■  6; 

n  =  1,2 . 6);  ^  " 

(b)  Symmetric  Longitudinal  Mode  (a  =  4068  psi;  k  =  15;  k  =  7; 

n  =  4,8,  .  28).  " 


Figure  4 


'J 

1  w'  4  ‘'^. 

different  from  one  another  but  the  buckling  stresses  are  almost  identical  - 
4068  psi).  This  buckling  stress  is  very  close  to  that  for  the  same  cylinder  without 
stiffeners  (cr  *  4047  psi;  ■  36613  pounds;  m  *  2;  n  »  3).  The  ’"esults  indicate 

that  small  number  of  discrete  stiffeners  is  not  effective  in  raising  the  buckling 
strength  {a^  *  4047  psi)  for  the  same  monocoque  cylinder  even  if  the  mode  shapes 
are  drastically  changed.  In  the  case  of  flat  plate,  however,  discrete  stiffening 
increases  the  buckling  strength  substantially. 

The  buckling  load  for  this  problem  was  also  calculated  based  on  the  uniform 
strain  assumption.  The  antisymmetric  circumferential  mode,  the  antisymmetric 
longitudinal  modes  (m  *  2,  4,  ...32),  and  the  circumferentially  related  modes  (n  * 

2,  6,  ..,.30)  governed  the  buckling  equations.  The  buckling  values  obtained  were: 

P  =  21632  pounds;  o  *  1639  psi;  a.  *  229.5  psi.  The  mode  shape  is  shown  in  Fig. 

X  X  o 

5.  It  is  of  interest  to  note  that  this  buckling  load  is  even  less  than  that 
(36613  pounds)  for  the  same  cylinder  without  stiffeners.  This  difference  is  caused 
by  the  basic  assumptions.  In  the  case  of  uniform  stress  assumption,  there  is  no 
circumferentical  stress.  But  in  the  case  of  uniform  strain  assumption,  the  ring 
introduces  high  compressive  circumferential  stress  in  the  skin  that  destabilizes 
the  shell . 

It  is  noted  that  the  computing  time  and  storage  were  drastically  reduced  in 
treating  the  submatrices  of  smaller  sizes  in  Eqs.  (28  and  29). 

3.  COMPARISON  WITH  EXPERIMENTAL  RESULTS 

Esslinger  and  Geier  (Reference  19)  tested  a  series  of  discretely  stiffened 
Mylar  cylinders  with  edges  cast  into  rigid  end  plates.  The  modulus  of  elas¬ 
ticity  of  Mylar,  which  is  not  available  in  Reference  19,  has  a  wide  range. 

However,  based  on  the  geometric  parameter  Z  of  5320,  the  Poisson's  ratio  of  1/3, 
and  the  theoretical  buckling  load  of  78  kg  given  in  Reference  19  for  a  monocoque 
cylinder,  the  modulus  of  elasticity  was  obtained  as  553.6  kg/mm^  by  using  the 
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classical  formula  that  the  Euler  buckling  stress  is  Eh/[3(l-v  )]  '  r.  Two  very 
stiff  rings  were  placed  at  the  ends  of  the  cylinder  to  simulate  the  built-in 
edge  conditions  that  w  =  v  =  3w/3x  =  =  0. 

(a)  Cylinder  With  Eight  Outer  Stringers  and  Two  Inner  Rings  Under  External 
Pressure 

The  sizes  of  the  cylinder  was  defined  as:  L=  330mm;  r=  100mm;  h=  0.193  mm. 

The  stiffeners  were  T-type  with  two  back-toback  angle  bars  of  6  x  6  x  0.193  mm. 

The  buckling  computations  were  based  on  the  antisymmetric  circumferential 
mode,  the  symmetric  longitudinal  modes  (m  =  1,  3,  ...41),  and  the  circumferentially 
related  modes  (n  =  4,  12,  20,  ...52).  For  the  uniform  strain  and  uniform  stress 
assumptions,  the  buckling  pressures  obtained  were  0.009768  and  0.01002  kg/cm^, 
respectively.  Both  are  in  close  agreement  with  the  experimental  value  of 
0.0095  kg/cm^.  The  buckling  mode  shape  was  plotted  in  Fig.  6  with  m  =  3  and 
n  =  12.  The  mode  shape  and  mode  numbers  agree  well  with  those  shown  in 
Fig.  11  of  Reference  19.  For  the  same  cylinder  without  stiffeners  the  buckling 
pressure  is  only  about  0.0039  kg/cm^  (Reference  19). 

(b)  Cylinder  With  Eight  Outer  Stringers  and  Two  Inner  Rings  Under  Uniform 
Axial  Compression 

The  cylindrical  shell  considered  was  the  same  as  the  previous  one  but  all 
the  stiffeners  were  single  angle  bars  of  6  x  6  x  0.193  mm. 

The  stiffeners  were  first  neglected.  In  the  buckling  computation,  the 
antisymmetric  circumferential  mode  with  n  =  6  and  the  symmetric  longitudinal 
modes  with  m  =  1,  3  ...79  were  used.  The  buckling  load  obtained  was  78.2  kg 
which  agreed  well  with  the  theoretical  value  of  78  kg  given  in  Reference  19. 

The  buckling  mode  shape  found  (with  one  half  longitudinal  sine  wave)  also  agreed 
with  the  experimental  photograph  (Reference  19). 


As  to  the  buckling  load  of  the  stiffened  shell.  Reference  19  provided  an 
experimental  value  of  45  kg,  a  theoretical  value  of  80  kg,  and  an  experimental 
postbuckling  value  greater  than  95  kg.  The  present  computations  were  based 
on  the  same  modes  as  those  used  in  the  case  of  external  pressure.  Based  on 
the  uniform  strain  assumption,  the  buckling  load  obtained  was  75.7  kg.  The 
mode  shape  is  shown  in  Fig.  7.  The  75.7  kg  predicted  is  close  to  the  theoretical 
value  of  30  kg  given  in  Reference  19.  The  present  buckling  mode  shape  agrees  well 
with  that  obtained  experimentally  in  Reference  19.  The  lack  of  agreement  between 
the  theoretical  and  experimental  buckling  loads  in  the  case  of  axially  compressed 
monocoque  cylindrical  shell  is  well-known.  This  phenomenon  also  appears  in  the 
present  case  of  discretely  stiffened  cylindrical  shell. 

A  buckling  load  of  92.0  kg  was  obtained  when  the  uniform  stress  assumption 
was  used.  This  value  should  be  expected  to  be  higher  than  that  (75.7  kg) 
obtained  by  using  the  uniform  strain  assumption.  The  reason  lies  mainly  in 
the  difference  in  circumferential  stress.  The  circumferential  stress  is  zero 
in  the  case  of  uniform  stress,  while  in  the  case  of  uniform  strain,  it  is  1.1% 
of  the  longitudinal  stress  which  corresponds  to  90. 8%  of  the  circumferential 
buckling  stress  for  the  same  cylinder  without  stiffeners  under  external  pressure. 

This  circumferential  stress  in  the  shell  is  introduced  by  the  two  rings  in 
circumferential  tension. 
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SECTION  V 


WEIGHT.  CONSTRAINTS,  AND  DESIGN  VARIABLES 

The  weight  W  of  an  orthogonally  stiffened  shell  with  multiple  stiffener 
sizes  is  expressed  as 

Weight  =  weight  of  the  skin  +  weight  of  the  smeared-out  stringers 
+  weight  of  the  smeared-out  rings  -  weight  of  the  material 
at  the  intersections  of  the  smeared-out  stringers  and 
rings  +  weight  of  the  discrete  stringers  +  weight  of  the 
discrete  rings  -  weight  of  the  material  at  the 
intersections  of  the  discrete  stringers  and  rings 
or 

W  =  2TrrLhp  +  27rrLhgjp  +  Tr[2r  -  sgn(d2)(d2+h)]hg2Lp 

-min(|dj  .[d^Dd  +  sgn(d^-d2)}/2-hg^hg2/|di*d2l  +  U  f 

k“  1  k 

+TT[2r  -  sgn(d  )(d  +h)]p  l\  -  min(  |d  I ,  |d  | ) 

k=l  '^k  ^ 

{1  +  sgn(d5-d^)}/2  tgt^n^n^p  (1) 

with 

h  =  d  t  /a 
ei  1  r  1 

sgn  (x)  =  x/|x| 

where  r,  h,  L  are  the  mean  radius,  thickness,  and  length  of  the  cylindrical 
shell,  respectively;  p  is  the  material  density,  d  and  t  are  the  depth  and 
width  of  the  stiffeners,  respectively;  the  subscripts  i,  2,  s,  and  r  are 
associated  with  the  smeared-out  stringer,  the  smeared-out  ring,  the  discrete 
stringer,  and  the  discrete  ring,  respectively;  a  is  the  spacing  between  the 
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stiffeners;  A  and  n  are  the  cross-sectional  area  and  the  total  number  of  the 
stiffeners  (See  Figure  1). 

For  the  cylindrical  shell  with  orthogonal  smeared-out  stiffeners  the 
inequality  constraints  are  defined  as  follows: 


(1)  buckling  of  the  stiffened  cylinder 


^  =  1  -  P  /p  >0 

x'  xo  - 


(42) 


where  is  the  general  buckling  load; 


(2)  buckling  of  the  stiffened  cylinder  between  the  circumferential 
rings  with  two  circular  edges  considered  as  simply-supported 

*2  '  ’  -  '’xp/'’xo2  («) 

(3)  buckling  of  the  shell  skin  bounded  by  longitudinal  and  circum¬ 
ferential  stiffeners  with  four  edges  considered  as  simply-supported 

Tr2h2E 


^3  = 


-  ^  0  -1  >  0; 
2  X  - 


3(l-v2)aj 

(4)  skin  failure  by  yielding 


(44) 


>  0 
y  - 


(45) 


where  is  the  yielding  stress  of  the  skin; 


(5)  longitudinal  stringer  failure  by  yielding 


1  -  :  0; 


(46) 


(6)  circumferential  ring  failure  by  yielding 


(47) 


(7)  buckling  of  the  longitudinal  stringer  between  the  circumferential 

rings  with  one  end  free  and  the  other  end  simply-supported( Reference  44 

.2_.d,  2 


'^7  "  170'- v2)  0.425]ra^j  -1  >  0; 


(48) 
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where  it  is  assumed  in  Eq.  (48)  that  both  stringers  and  rings  are 
on  the  same  side  of  the  skin  and  that  the  rings  are  deeper  than  the 
stringers.  If  these  conditions  are  violated,  the  spacing  ^2 
Eq.  (48)  should  be  replaced  by  the  spacing  between  two  discrete  rings 
L/(l  +  n^)  on  the  same  side  of  the  skin.  In  the  absence  of  discrete 
rings,  n^  is  zero; 

(8)  design  variable  constraints 

<  u.  <  u  i  =  1 ,2, . . .  p  (49) 

i.  -  ^  -  u^ 

O'"  .  , 

t|;  =  (u.  -  )(u^  -  u.)>  0  (50) 

where  the  design  variables  consist  of  a^,  a^,  t^,  t^, 
dj^,  d^,  and  h. 

The  design  of  the  cylindrical  shell  with  orthogonal  smeared-out 
stiffeners  thus  includes  seven  design  variables  and  fourteen  inequality 
constraints. 

For  the  cylindrical  shell  with  orthogonal  smeared-out  and  discrete 
stiffeners,  the  inequality  constraints  described  for  the  smeared-out 
case  still  hold  but  the  computation  of  the  allowable  uniform  axial  buck¬ 
ling  load  used  in  Eqs.  (3)  and  (4)  should  incorporate  the  effect  of 
discrete  stiffeners  as  well.  Besides,  three  more  constraints  are  included: 

(1)  longitudinal  discrete  stringer  failure  by  yielding 

(2)  circumferential  discrete  ring  failure  by  yielding 

'I'  =  1  -  l^erl^V  -  °  ^ 
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(3)  buckling  of  the  longitudinal  discrete  stringer  between  two 
circumferential  rings  with  one  end  free  and  the  other  end 


simply  supported 
t. 


E 


"I- 


=  12(1-v2T  *3:’  '■*  1- 


d-(l+n„)  2 


■}  +  0.425] 


a  -  1  >  0 

xs 


(53) 


It  is  assumed  in  Eq.  (53)  that  both  stringers  and  rings  are  on  the 
same  side  of  the  skin  and  that  the  rings  are  deeper  than  the  stringers. 

If  these  conditions  are  violated,  n^  should  be  set  equal  to  zero. 

The  design  variables  chosen  for  this  case  are  a  ,  a  ,  t  ,  t„,  d  , 
d^,  h,  tg,  dg,  t^  and  d^.  The  design  of  the  cylindrical  shell  with  both 
smeared-out  and  discrete  stiffeners  thus  includes  11  design  variables  and 
21  inequality  constraints. 
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SECTION  VI 

A  SIMPLIFIED  BUCKLING  FORMULATION 

It  has  been  shown  in  the  example  2  of  Section  IV  that  for  the 
buckling  analysis  of  a  simply  supported  cylindrical  shell  with  one  ring 
and  four  stringers  under  uniform  axial  compression,  the  exact  formulation 
required  a  minimum  of  750  eigenvalue  equations  (with  =  25  and  = 

k^  =  5) .  With  the  use  of  a  i.iatrix  uncoupling  technique  based  on  the  as¬ 
sociated  mode  shapes,  a  matrix  reduction  technique  by  partitioning  and 
substitution,  a  special  compact  storage  scheme  to  make  use  of  the  sparse¬ 
ness  of  the  submatrices,  and  the  Ritz  iteration  method  combined  with 
the  Chebyshev  eigenvalue  solution  procedure,  the  lowest  eigenvalue  for  the 
750  equations  was  able  to  be  solved  eventually  by  a  CDC  6500  computer 
with  only  approximately  400  seconds  of  central  processing  time.  Although 
it  has  already  been  a  tremendous  reduction  in  computation,  400  seconds  are 
still  formidably  long  if  the  buckling  load  has  to  be  computed  repeatedly 
for  optimization. 

To  circumvent  this  difficulty,  the  computing  time  for  buckling  load 
must  be  further  reduced  drastically  and  it  appears  that  one  has  no  al¬ 
ternative  but  to  seek  for  an  extremely  simple  and  short  method  for  only 
rough  estimate  of  buckling  load.  In  that  case  the  exact  buckling  compu¬ 
tation  still  has,  of  course,  to  be  performed  for  the  optimized  design 
variables  to  check  the  validity  of  the  buckling  load  obtained  by  the 
approximate  method. 


r 

( 
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In  the  case  of  cylinder  with  smeared-out  stiffeners,  it  was  shown  in 
Reference  10  that  an  exact  solution  can  be  achieved  by  using  one  term  for 
each  u,  v,  and  w  functions.  Thus  for  extreme  simplicity,  one  anti¬ 
symmetric  circumferential  mode  and  one  longitudinal  mode  are  assumed  in 
the  displacement  functions  for  the  cylinder  with  discrete  stiffeners. 

f  u  =  cos  ~-sin  ne 
mn  L 

J  V  =  V„  sin  cos  ne  (54) 

\  mn  L  '  ' 

w  =  W  sin  sin  ne 
1.  mn  L 

The  displacement  functions  are  substituted  into  the  potential 
energy  expressions  (Eq.  (1))  and  the  same  procedure  as  shown  in 
Section  II  is  followed  for  deriving  the  buckling  equations.  The  resulting 
equations  are  similar  to  the  ones  for  the  cylinder  with  smeared-out  stif¬ 
feners  and  its  explicit  formulations  can  be  obtained  by  setting  m=m^.=mj 
and  n=n^=nj  in  Eq.(27)  .  The  critical  buckling  load  is  obtained 
by  searching  for  the  smallest  eigenvalue  in  a  set  of  three  equations 
by  varying  the  values  of  integers  m  and  n. 

In  the  case  of  cylinders  with  smeared-out  stiffeners  a  smooth  sur¬ 
face  for  the  minimum  eigenvalue  is  obtained  so  that  the  efficient  section¬ 
ing  approach  proposed  and  demonstrated  in  Reference  37  may  be  used  for  the 
search  for  the  minimum.  In  the  case  of  cylinders  with  discrete  stiffeners, 
there  are  many  local  minima  in  the  eigenvalue  surface.  The  absolute  min¬ 
imum  has  to  be  searched  for  by  varying  m  and  n  continuously  in  a  certain 
range.  In  spite  of  this  disadvantage,  this  simple  approximate  discrete 
calculation  is  superior  to  the  smeared-out  calculation  in  the  buckling 
prediction  of  cylinder  with  discrete  stiffeners. 


This  approximate  discrete  method  was  first  verified  by  performing 
a  buckling  computation  of  a  cylinder  with  100  stringers  and  16  rings. 

The  material  and  geometry  parameters  of  the  cylinder  are  defined  as: 

E  =  l.OSxlO^psi;  v  =  1/3;  a-)  =  a2  =  3  in.;  d]  =  d2  =  0.2  in.;  t]  =  t2  = 
0.125  in.;  h  =  0.05  in.;  r  =  48  in.;  and  L  =  50  in.  The  exact  solution 
obtained  by  the  exact  discrete  formulation  was  given  in  Section  IV.  The 
buckling  load  was  801.3  Ibs/in  and  the  buckling  mode  numbers  were  m  =  6 
and  n  =  17.  The  present  approximate  discrete  computation  gave  a  buck¬ 
ling  load  of  807.9  Ibs/in  with  the  mode  numbers  m  =  6  and  n  =  17.  It  is 
noted  that  in  order  to  obtain  the  converged  solutions  by  using  the  dis¬ 
crete  formulations,  600  degrees  of  freedom  were  required  while  in  this 
approximate  method  only  three  degrees  of  freedom  were  used. 

This  approximate  discrete  method  was  then  tested  for  a  cylinder 

with  only  one  ring  and  four  stringers.  The  material  and  geometry  parameters 

of  the  cylinder  are  defined  as:  E  =  10.6x10®  psi ;  v  =  1/3;  h  =  0.03  in.; 

r  «  48  in.;  L  =  50  in.;  d  =  d  =  1.3196  in.;  t  =  t„  =  0.8248  in.  The 

s  r  s  I 

exact  converged  buckling  stress  was  given  in  Section  IV  as  4068  psi  with 
several  different  mode  shapes.  The  present  discrete  approximate  method 
gave  a  buckling  stress  of  5008  psi  with  buckling  mode  numbers  m  =  24  and 
n  =  2. 

It  is  important  to  note  that  when  the  seven  sets  of  near  minimum 
eigenvalue  solutions  fo>  the  above  examples  were  obtained  in  Section  IV, 
the  sizes  of  the  originally  coupled  eigenvalue  equations  were:  504;  630; 

750;  936;  1008;  1134;  and  1248.  If  the  present  approximate  method  is 
used,  only  three  eigenvalue  equations  need  be  solved.  The  loss  in  ac¬ 
curacy  is  compensated  for  by  the  drastic  reduction  in  computing  time. 

Without  such  drastic  reduction  the  optimization  computations  will  be  too 
expensive  to  perform.  The  buckling  load  for  the  optimized  design 
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variables  obtained  by  this  approximate  method  is,  however,  checked  by 
the  expensive  exact  discrete  computation.  The  difference  in  the  pre¬ 
dicted  buckling  loads  between  this  approximate  discrete  method  and  the 
exact  discrete  method  becomes  smaller  as  the  number  of  stiffeners  is 
increased.  As  will  be  shown  in  a  subsequent  example  of  minimum  weight 
design,  the  buckling  load  predicted  by  this  approximate  method  for  a 
cylinder  with  5  discrete  rings  and  30  discrete  stringers  is  very  close 
to  the  exact  discrete  solution. 


SECTION  VII 


THEORY  OF  OPTIMIZATION 

The  theory  of  optimization  has  long  been  an  attractive  research 
topic  and  extensive  research  works  have  been  done.  Various  theories 
of  optimization  have  been  applied  to  the  minimum  weight  design  of  struc¬ 
tures  with  various  constraints.  These  applications  have  been  discussed 
in  the  Introduction. 

In  Reference  40,  Schuldt  developed  a  method  of  multipliers  for  mathe¬ 
matical  programming  problems  with  equality  and  inequality  constraints. 

Since  this  method  has  not  been  applied  to  the  minimum  weight  design  of 
structures,  it  is  introduced  here  as  an  efficient  method  for  the  min¬ 
imum  weight  design  of  the  stiffened  cylindrical  shell  with  multiple 
stiffener  sizes. 

The  conventional  problem  of  optimum  design  is  to  minimize  the 
weight  of  a  structure  expressed  as  a  function  of  the  design  variables 
which  satisfy  a  series  of  constraints  on  the  design  behavior  and  on  the 
sizes  of  the  variables  (design  variable  equality  constraints  or  design 

variable  inequality  constraints).  That  is,  to  minimize 

W(u)  (55) 

subject  to  ^(u)  =  0 

3;(u)=0  (56) 

According  to  Reference  40  the  problem  may  be  transformed  into  the  one 
of  minimizing  the  augmented  function 

=  U(u)  +  f  W)  -  (57) 
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subject  to 

^(u)  =0  (58 ) 

^(u)  =0  (59) 


which  satisfy  the  following  first-order  condition  at  the  solution  point. 


Fy(u,X,y)  =  W^j(u)  +  ({>y(u)X  -  <j;jj(u)y  =  0 


V*,,.  'J^.,.(u)  0  j  l»»**»r 

(3)  (:) 


j— 1 t...yr 


(60) 

(61) 

(62) 


The  augmented  penalty  function  is  written  as 


^(uj.y.k)  =  Q(u;^.k)  +  R(i(,?;.k) 

where 


(63) 


Q  =  W(u)  +  f?(u)  +  k  |^(u)^(u) 

=  W(ii)  +  k[?(2f)  +  ^  [^(^()  +  ^  5:]  -  ^  ft  (64  ) 

R  =  k[j(u)  -  ^  u]^[^(u)  -  ^  v*]  -  ft  (65 ) 

i|;j(u)  =  minU^jj(u),  ^  M^jj]  j=l,....r  (66) 

The  constant  k  is  initially  chosen  from  scaling  consideration  and  held 
fixed  throughout  the  algorithm. 

r\j  f\j 

Initial  cycle  is  begun  with  the  choice  of  X  =  0  and  v  =  0,  i.e., 
fi(u.0,0,k)  =  W(u)  +  k[f(u)^(u)  +  f(u)3;(u)]  (67) 

where 

illj(x)  =  min{i;;j(x),0}  j=l....,r  (68) 
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The  following  multiplier  updating  rules  are  applied. 


5:,  =  5r,  1 2k  i(x)  (69) 

=  S,  -  2k  Kx)  (20) 

where  subscripts  1  and  2  correspond  to  the  current  and  the  next  cycles, 
respectively. 

Differentiating  Eq.  (63)  with  respect  to  the  design  variables 
u,  the  following  necessary  conditions  for  the  minimum  is  obtained. 


(u ,A,vi,k)  - 


w,(u)  +  +  2k  Wn  -  *(u)[y  -  2k  =  0  (71) 


Combining  Eqs.  (  69),  (  70),  and  (71),  one  may  obtain  the  following 
expression  where  the  updated  multipliers  satisfy  the  optimality  con¬ 
dition  given  by  Eq.  ( 60) 


yj  =  w. 


(72) 


In  addition  to  this  condition,  the  equality  constraints  by  Eq. 
(58)  and  the  inequality  constraints  defined  by  Eq.  (59)  should  be 
satisfied.  In  Reference  40  Schuldt  first  showed  that  the  augmented 
penalty  function  increases  at  the  end  of  each  minimization  cycle  when 
the  multipliers  are  changed  from  and  to  A^  and  respectively. 
He  then  demonstrated  through  one  example  that  succeeding  cycles  of 
minimizing  the  augmented  penalty  function  may  drive  toward  constrained 
satisfaction. 

In  order  to  minimize  the  augmented  penalty  function  defined  by 
Eq.  (63),  the  method  for  the  unconstrained  function  minimization  oro- 
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posed  by  Fletcher  and  Powell  (Reference  45)  and  a  quasi -Newton  method 
are  used.  The  two  methods  are  also  compared. 

All  gradient  calculations  are  carried  out  by  the  simple  finite 
difference  method  in  this  study. 


SECTION  VIII 


TWO  OPTIMUM  WEIGHT  DESIGN  EXAMPLES 

The  first  design  example  chosen  was  a  simply-supported  orthogonally 
stiffened  cylindrical  shell  subjected  to  uniform  axial  compression.  The 
cylinder  had  a  radius  of  48  inches  and  a  length  of  50  inches.  The  pro¬ 
perties  of  the  aluminum  were  used  for  both  the  cylinder  and  the  stiffeners: 
E  =  10.6xl0®psi;  v  =  1/3;  p  =  0.101  Ibs/in^.  The  buckling  load  of 
P^Q  =  238.2  kips  or  =  789.8  Ibs/in  was  used  as  the  buckling  constraint. 
The  seven  initial  design  variables  which  were  restrained  between  0.01 
and  10  inches  are  shown  in  Table  4.  This  problem  is  similar  to  the  one 
given  in  Reference  37  but  with  one  extra  design  variable  ^2' 

The  buckling  computation  was  based  on  the  simple  smeared-out  formu¬ 
lation  and  the  optimization  process  was  performed  by  the  method  of  mul¬ 
tipliers  with  the  Fletcher-Powell  method  as  well  as  the  quasi-Niwton 
method.  One  of  the  purposes  of  choosing  this  example  was  to  obtain 
data  for  comparison  with  a  later  example  with  additional  discrete  stif¬ 
feners.  The  other  purpose  of  conducting  this  example  was  to  evaluate 
the  performance  of  the  present  optimization  methods.  The  results  for 
the  minimum  weight  design  are  presented  in  Table  4.  For  this  example  the 
Fletcher-Powell  method  seems  to  be  more  suitable  than  the  quasi -Newton 
method. 

The  second  example  was  a  simply-supported  aluminum  cylindrical  shell 
with  5  discrete  rings,  30  discrete  stringers,  and  smeared-out  stiffeners. 
The  cylinder  had  the  same  length,  radius,  and  material  properties  as  the 
first  example.  The  buckling  constraint  was  789.82  Ibs/in.  The  11  initial 
design  variables  which  were  restrained  between  0.01  and  10  inches  are 
shown  in  Table  5. 
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Table  4  Weight  minitnization  using  the  smeared-out  formulation 

(Nxo  =  789.82  Ibs/in.;  =  50  ksi;  0.01  in.  <  <  10  in., 

i=l,2,...7) 


Mode 

Number 


913.5 

789.5 


4  15 
6 


3.0 

4.962 


Method 


.2 

■E 

.2 

0.125 

0.125 

0.05 

0 

Fletcher 

.273 

li 

.7 

0.02287 

0.00999 

0.01559 

2315 

Powel 1 

0.125 

0.02546 


0.125 

0.01181 


0.05 

0.02126 


quasi - 
Newton 


Table  5  Weight  minimization  using  smeared-out  and  discrete  formulations 
(N..»  =  789.82  Ibs/in.;  a,  =  50  ksi;  0.01  in.  <  u.  <  10  in.  with 


Iteration  Weight  N  Mode  a,  a  d 

NFE*  ^  T  ^  ^ 

er  (lbs.)  (Ibs/in)  m  n  (in.)  (in.)  (in.) 


0 

141.21 

1102.8 

3 

1636 

42.83 

815.8 

11 

0.2000 

1.0446 


0  141.21 

34  78.01 


0  3.0000  0.2 
9  3.0003 


t,  h 

(in.)  (in.)  (in.) 


t  d  t 

s  r  r 

(in.)  (in.)  (in.) 


Method 


0.1250 

0.0500 

0.4085 

zzrzj 

0.2553 

0.0913 

0.0110 

0.0100 

0.6938 

0.0896 

0.0145 

0  Fletcher- 
1711  Powell 


0.0500 

0.0230 


0.2553 

0.2741 


0.0571  0  quasi- 

0.0509  841  Newton 


*  Number  of  function  evaluations 
t  CDC  6500  central  processing  time 


! 

t 

The  buckling  computation  during  the  weight  minimization  was  per¬ 
formed  using  the  present  approximate  discrete  technique.  The  results 
of  the  minimum  weight  design  are  given  in  Table  5.  The  approximate 
buckling  load  for  the  minimum  weight  design  obtained  by  using  the 
Fletcher-Powell  method  was  also  checked  by  using  the  exact  buckling 
computation  based  on  the  discrete  formulation.  The  symmetrical  longi¬ 
tudinal  modes  with  m  =  1,  3,  5,  ...  59,  the  antisymmetric  circumferential 
modes,  and  the  circumferentially  related  modes  with  n  =  7,  23,  37,  53 
governed  the  buckling  equations.  The  exact  buckling  load  obtained  was 
801.7  Ibs/in  which  was  very  close  to  the  815.8  Ibs/in  obtained  by  the 
approximate  discrete  method.  The  corresponding  mode  shape  is  shown  in 
Fig.  8.  This  optimized  mode  shape  shows  that  a  very  smooth  buckling  surface 
resulted  from  the  combined  effect  of  the  smeared-out  stiffeners  and  the 
discrete  stiffeners. 

Again,  Table  5  reveals  that  the  Fletcher-Powell  method  is  more 
suitable  in  the  present  weight  minimization  than  the  quasi-Newton  method. 

The  process  of  the  weight  minimization  by  using  the  Fletcher-Powell 
method  is  summarized  in  Table  6.  The  weight  reduction  was  completed  in 
the  early  stage  of  the  optimization  but  the  violation  of  the  constraints 
was  sometimes  greater  than  ten  percent.  The  rest  of  the  optimization 
cycles  was  carried  out  to  correct  the  violations.  All  the 
final  designs  violated  the  constraints  by  less  than  one  percent. 
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The  optimized  stiffened  shell  with  single  stiffener  size  has  the 
failure  modes  of  general  buckling  and  buckling  of  the  shell  skin  bounded 
by  longitudinal  and  circumferential  stiffeners,  expressed  by  Eq.  (42)  and 
Eq.  (44),  respectively.  The  optimized  stiffened  shell  with  multiple 
stiffener  sizes  is  inside  the  feasible  region  and  near  the  failure  modes 
of  general  buckling,  buckling  of  the  shell  skin,  and  buckling  of  the 
smeared-out  longitudinal  stringer,  expressed  by  Eq.  (42),  Eq.  (44)  and 
Eq.  (48),  respectively. 

The  present  examples  show  that  the  method  of  multipliers  with  the 
Fletcher-Powell  method  is  very  effective  in  the  weight  minimization. 

The  results  in  Table  6  show  that  under  the  same  buckling  constraint  the 
cylinder  with  smeared-out  and  discrete  stiffeners  results  in  less  weight 
than  the  ordinary  cylinder  with  only  smeared-out  stiffeners. 
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SECTION  IX 


CONCLUDING  REMARKS 

The  buckling  equations  for  the  cylindrical  shell  with  multiple  stif¬ 
fener  sizes  are  formulated  using  the  energy  method.  For  the  case  of  equally 
spaced  and  identical  stringers  and  equally  spaced  and  identical  rings,  a 
close  examination  of  the  buckling  equations  leads  to  the  uncoupling  of 
the  buckling  equations  into  separate  sets  of  equations  associated  with 
the  symmetric  and  antisymmetric  longitudinal  modes,  symmetric  and  anti¬ 
symmetric  circumferential  modes,  and  circumferentially  related  modes.  By 
using  such  uncoupling  and  some  matrix  partitioning  and  reduction,  the  huge 
sets  of  matrix  equations  required  for  the  buckling  predictions  of  a  com¬ 
mon  discretely  stiffened  cylinder  can  be  reduced  drastically  to  manageable 
sizes.  Effort  is  made  to  preserve  the  sparseness  of  the  buckling  matrices 
so  that  a  special  compact  storage  technique  can  be  used. 

A  method,  the  Ritz  iteration  with  Chebyshev  procedure,  is  developed 
for  efficient  computation  of  the  minimum  eigenvalue  of  a  large  general 
eigenvalue  problem.  The  superiority  of  this  method  over  other  methods  is 
demonstrated  through  a  numerical  example. 

The  discrete  formulations  are  first  used  to  analyze  an  axially  compressed 
cylinder  with  a  large  number  of  stiffeners  for  which  an  alternative  solution 
based  on  the  smeared-out  method  can  be  obtained  for  comparison.  The 
agreements  are  very  good  for  the  formulations  based  on  both  the  uniform 
stress  and  uniform  strain  assumptions  (see  Table  2).  The  formulations 
are  then  verified  by  comparing  the  computed  buckling  loads  and  modes  with 
the  experimental  results  for  a  clamped  cylinder  with  two  rings  and  eight 
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stringers  under  external  pressure  and  then  under  axial  compression.  For 
the  case  of  external  pressure,  the  agreements  are  good  for  the  formulations 
based  on  both  the  uniform  stress  and  uniform  strain  assumptions.  For 
the  case  of  axial  compression,  the  agreements  are  poor.  Such  lack  of 
agreement  between  theoretical  and  experimental  axial  buckling  loads  has 
long  existed  in  the  case  of  monocoque  cylinder. 

It  should  be  noted  that  both  the  uniform  stress  and  uniform  strain 
assumptions  are  shown  to  be  consistent  with  each  other  and  correct  for 
the  cylinder  under  external  pressure.  This  is  the  important  case  where 
a  small  number  of  discrete  stiffeners  is  found  to  be  very  effective  in 
raising  the  buckling  loads.  Such  conclusion  was  made  previously  by 
Esslinger  and  Geier  in  an  experimental  study.  For  the  case  of  axial 
compression,  the  two  prestress  assumptions  do  not  result  in  the  same 
buckling  load.  The  uniform  strain  assumption  results  in  circumferential 
compression  in  the  skin  and  consequently,  results  in  lower  buckling 
load  than  the  uniform  stress  assumption.  However,  this  is  the  unimportant 
case  where  a  small  number  of  discrete  stiffeners  is  found  to  be  ineffective 
in  raising  the  buckling  loads. 

The  minimum  weight  design  of  orthogonally  stiffened  cylindrical 
shells  with  discrete  stiffeners  in  addition  to  the  smeared-out  stiffeners 
has  not  been  studied  before  mainly  due  to  the  computational  difficulty 
in  obtaining  the  buckling  load.  With  the  development  in  Sections  II,  III, 
and  IV,  which  formulated  the  discrete  stiffeners  exactly  and  simplified 
the  computational  method  tremendously,  such  computational  difficulty  was 
reduced.  The  simplifications  included  uncoupling  the  buckling  matrix  based 
on  the  associated  mode  shapes,  reducing  the  uncoupled  submatrices  by 
partitioning  and  substitution,  taking  advantage  of  the  sparseness  of  the 
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buckling  matrix  by  a  special  compact  storage  scheme,  and  using  an  ef¬ 
ficient  eigenvalue  solution  procedure.  The  development  Is,  however, 
still  too  expensive  to  be  used  repeatedly  for  weight  minimization. 

In  this  report,  an  approximate  discrete  formulation  based  on  very 
simple  displacement  functions  has  been  suggested  for  computing  the  buck¬ 
ling  load  for  the  discretely  stiffened  cylinders.  This  method  requires 
the  solution  of  a  set  of  only  three  Instead  of  several  hundred  eigenvalue 
equations.  The  difficulty  In  computing  time  for  the  weight  minimization 
has  thus  been  overcome.  To  ensure  that  the  optimized  design  using  this 
approximate  discrete  method  for  predicting  buckling  loads  Is  valid,  the 
exact  discrete  buckling  computation  Is  still  performed  based  on  the 
final  design  variables. 

To  cope  with  this  simple  discrete  method  for  buckling  computation, 
the  methods  of  multipliers  with  the  Fletcher-Powell  method  and  with 
quasi-Newton  method  have  been  chosen  for  the  weight  optimization.  Two 
design  examples  have  been  performed. 

The  design  examples  have  provided  the  following  Indications: 

(1)  The  present  simple  discrete  method  Is  very  effective  In  weight 
minimization  design.  The  second  example  shows  that  the  buckling  load 
predicted  by  this  approximate  method  was  very  close  to  the  exact  solu¬ 
tion. 

(2)  The  comparison  shows  that  for  unconstrained  minimization  the 
Fletcher-Powell  method  Is  more  suitable  than  the  quasi -Newton  method. 

The  former  method  Is  very  effective  In  dealing  with  the  present  type 
of  weight  minimization  problem. 
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(3)  For  the  minimum  weight  design  of  cylindrical  shell  under  buck¬ 
ling  constraints,  the  cylinder  with  both  the  smeared-out  and  the  discrete 
stiffeners  can  be  lighter  than  the  cylinder  with  only  smeared-out 
stiffeners. 

Finally,  it  should  be  noted  that  the  above  indications  are  obtained 
from  only  two  design  examples  and  they  should  not  be  regarded  as  g.ineral 
conclusions . 
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